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Abstract 

For some centuries, first order chemical rate constants were determined mainly 
by a linear logarithmic plot of reagent concentration terms against time where the 
initial concentration was required, which is experimentally often a challenging task 
to derive accurate estimates. By definition, the rate constant was deemed to be in- 
variant and the kinetic equations were developed with this assumption. A reason 
for these developments was the ease in which linear graphs could be plotted. Here, 
different methods are discussed that does not require exact knowledge of initial 
concentrations and which require elementary nonlinear analysis and the ensuing 
results are compared with those derived from the standard methodology from an 
actual chemical reaction, with its experimental determination of the initial concen- 
tration with a degree of uncertain. We verify experimentally our previous theoret- 
ical conclusion based on simulation data [ J. Math . Chem 43 (2008) 976-1023] 
that the so called rate constant is never constant even for elementary reactions and 
that all the rate laws and experimental determinations to date are actually aver- 
aged quantities over the reaction pathway. We conclude that nonlinear methods in 
conjunction with experiments could in the future play a crucial role in extracting 
information of various kinetic parameters. 

Keywords: [1] elementary reaction rate constant , [2] activity and reactivity coeffi- 
cients, [3] elementary and ionic reactions without pre-equilibrium. 
AMS Classification: 80A10, 80A30, 81T80, 82B05, 92C45, 92E10, 92E20. 

1 1. INTRODUCTION AND METHODS 

As alluded in the abstract, most kinetic determinations use logarithmic plots with 
known initial concentrations, although there have been attempts (fl~|[2] and refs. therein] . 
(There are possible ambiguities in [T) concerning choice of variables that will not be 
discussed.) However all these publications hitherto assume constancy of the rate con- 
stant k and do not focus on nonlinear analysis (NLA), as will be attempted here in 
preliminary form. We analyze kinetic data of the tert butyl chloride hydrolysis reac- 
tion in ethanol solvent (80%v/v) derived from the Year III teaching laboratory of this 
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University (UM); 0.3mL of the reactant was dissolved in 50mL of ethanol initially. 
The reaction was conducted at 30°C and monitored over time (minutes) by measuring 
conductivity (pS cm -1 ) due to the release of H + and CP ions as shown below (Q]), 

C 4 H 9 C1 + H 2 C 4 H 9 OH + H + +Cl (1) 

and Aqo = 2050^S cm -1 was determined by heating the reaction vessel at the end of 
the monitoring to 60°C until there was no apparent change in the conductivity when 
equilibrated back at 30°C. "Units" in the figures and text refers to /iS cm~ 1 . It would be 
inferred here that either because of evaporation or the temperatures not equilibrating 
after heating, the measured Aoo is larger than the actual one. Linear proportionality 
is assumed in A and the extent of reaction x, where the first order law (c being the 
instantaneous concentration and a the initial concentration) is 4| = — kc — —k(a — x); 
with Aoo = oca, X t = ax and A(0) = A = ax , integration yields for assumed 
constant k 

Eqn.(f2| determines k if Aq and A^, are known. The plot of (f2J) was made for the same 
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Figure 1 : Integrated equation© plot with from experiment (0) and from two dif- 
ferent arbitrary values (1,2) for Aoo 

experimental values with different Aoo 's, both higher and lower than the experimental 
value. We find that the rate constant for the NLA was higher, leading to a lower value 
of Aoo which is consonant with evaporation of solvent or the non-equilibration of tem- 
perature prior to measurement to determine Aoo -Except for the last subsection, we shall 
do a NLA based on constant k assumption. 



1.1 1.1 Method 1 

Under linearity argument and constant k, the rate equation 4| = — kc = —k(a — x) 
reduces to 

M = _fcA(t)+A 0o .fe (3) 

Hence a plot of vs t would be linear. We find this to be the case for polynomial 
order npoly < 3 as in Fig.© below for all data values; higher polynomial orders can 
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be used in selected data points of the curve below, especially in the central region. Thus 
criteria must be set up to determine the appropriate regime of datapoints in the NLA. 
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Figure 2: Method 1 graph showing linearity lower order polynomial fits 



1.2 1.2 Method 2 

Let a' = Aqo — An, then In a' — ln(Aoo — A) = kt, then noting this and differentiating 
yields 

In = -kt + ln[fc(Aoo - A )] (4) 

Y C 

A typical plot that can extract k as a linear plot of ln(d\/dt) vs t is given in Fig.©. 
Linearity is observed for npoly — 2 and smooth curves without oscillations for at least 

npoly < 3. 
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Figure 3: Method 2 where smooth curves are obtained for at least npoly < 4 
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1.3 1.3 Other methods and considerations 

A variant method similar to the Guggenheim method 12) of elimination is given below 
but where gradients to the conductivity curve is required, and where the average over 
all pairs is required. 

N N 

{k) = N(N -i) T,Y, ln ( A/ ^)/ A/ fe)) /(** - *i) (5) 

' i j>i 

It was discovered that the normal least squares polynomial method using Gaussian 
elimination [4, Sec.6.2.4,p.318 ]to derive the coefficients of the polynomial was highly 
unstable for npoly > 4 and so for this work, we used a variant of the Orthogonal 
method modified for determination of differentials. The normal method defines the 
nth order polynomial p n (t) which is then expressed as a sum of square terms over the 
domain of measurement to yield Q in eqns©. 

pn(t) = ej ,><r 2 (6) 

Q(f,Pn) = J2i^l[fi ~ Pn(U)] 

The Q function is minimized over the polynomial coefficient space. In the Orthog- 
onal method adopted here, we express our polynomial expression p m (i) linearly in 
coefficients dj of ipj functions that are orthogonal with respect to an inner product def- 
inition. For arbitrary functions /, g, the inner product (/, g) is defined below, together 
with properties of the ipj orthogonal polynomials. 



(/,<?) = ££=i/(**).fl(**) 

ipi(t) = (t - &j)</j,_i(i) - Ci(pi- 2 (t) (i > 1) 

<p (t) = l,and (fj = j < I, 

hi = (t</3j_i,y)j_i)/(^i_i,^i_i) (i > 1) 

Ci = (ttpi-i,(pi- 2 )/(<Pi-2, <pi-s) (i > 2),c< = 0. 

We define the m th order polynomial and associated eij coefficients as: 



(7) 



(8) 



(9) 



a 3 = {f,<Pj)/(<Pj><Pj)i(J = Q> 1 >--- m ) 

The recursive definitions for the first and second derivatives are given respectively 

as: 

¥>*(*) = <d-iW{t-bi)+<Pi-i(t)-Ci<d-2 (*)(<>!) rim 
v »{t) = ^ / _ 1 (t)(t-6 i ) + 2^_ 1 (t)-c i ^'_ 2 (*)(i>2) 

Here the codes were developed in C/C++ which provides for recursive functions which 
we exploited for the evaluation of all the terms. The experimental data were fitted to 
an m th order expression A m (t) defined below 

n 

Figure(|4]iare plots for the different polynomial orders n. The orthogonal polynomial 
method is stable and the mean square error decreases with higher polynomial order 
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Figure 4: Plot using orthogonal polynomials for various orders n 



(for the 36 data points) monotonically, but the differentials are not so stable, as shown 
in the previous figures. 

Differentiating © for constant k leads to ( TTZi i expressed in two ways 



fx 

dt 2 



= -k 



(a) or k 



d 2 X fdX 
"dl 2 ' \~dt 



(b) 



(12) 



Eq.(Q2tb)) suggests another way of computing k for "well-behaved" values of the 
differentials, meaning regions where k would appear to be a reasonable constant. The 
(a) form suggests an exponential solution. Define 4| = dl and = d2l. Then 
dl(t) = Aexp(-kt) and dl(Q) = A = h,2 from (flTT i . Furthermore . as t — * 0, k = 
{—Ihifhi) and a global definition of the rate constant becomes possible based on the 
total system X(t) curve. 

With a slight change of notation, we now define dl and d2l as referring to the 
continuous functions dl(t) = Acxp(-kt) and d2l(t) = —kAex-p(-kt) and we con- 
sider (dX/dt) and d 2 X/dt 2 to belong to the values (fTTt derived from Is fitting where 
(dX/dt) = X' m , (d 2 X/dt 2 ) = A"j which are the experimental values for a curve fit 
of order m. From the experimentally derived gradients and differentials, we can define 
two non-negative functions R a (k) and Rb(k) as below: 



Ra(k) = 

Rb(k) = 

where 



Ef =1 ( 



d\(tj) 
dt 



dl(U) 



(13) 



f a (k) = R' a (k) and f b (k)=R' b (k) 



and a minimum exists at f a (k) = fb(k) = 0. We solve the equations f a , /{, for their 
roots in k using the Newton-Raphson method and compute the rate constant k. The 
error threshold in the Newton-Raphson method was set at e = 1.0 x 10 ~ 7 We provide 
a series of data of the form [n, A, k a , kt,, X a .oo, X by0 o] where n refers to the polynomial 
degree, A the initial value constant as above, k a and k b is the rate constant for function 
f a and ft, (solved when the functions are zero respectively ) and likewise for A Q ,oo 
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and Xb,oo- The e symbol refers to base 10 (decimal) exponents. The A^ values are 
averaged over all the (36) data points from the equation 

A» = ^i + A(t) (14) 
at k 

The results are as follows: 

[2, 3.7634e0, 3.2876e - 3, 3.2967e - 3, 1.1506e3, 1.1477e3], 
[3, 3.6745e0, 2.7537e - 3, 2.7849e - 3, 1.34756e3, 1.3334e3], 
[4, 3.6380e0, 2.0973e - 3, 2.4716e - 3, 1.7408e3, 1.4900e3], 
[5, 4.0210e0, 9.7622e - 3, 4.9932e - 3, 4.4709e2, 7.9328e2, ], 
[6, 4.5260e0, 4.1270e - 2, 8.9257e - 3, 1.7101e2, 4.8403e2]. 

We noticed as in the previous cases that the most linear values occur for 1 < n < 4. 
In this approach, we can use the f a and /& function similarity of solution for k to 
determine the appropriate regime for a reasonable solution. Here, we notice a sudden 
departure of similar value between k a and fcf, (about 0.4 difference ) at n = 4 and so 
we conclude that the probable "rate constant" is about the range given by the values 
spanning n = 2 and n = 3. Interestingly, the values are approximately similar 
to the ones for method 1 and 2 for polynomial evaluation 2 and 3 for those methods. 
More study with reliable data needs to be done in order to discern and select appropriate 
criteria that can be applied to these non-linear methods. 

1.4 1.4 Evidence of varying kinetic coefficient k 

Finally, what of direct methods that do not assume the constancy of k which was the 
case in the above subsections? Under the linearity assumption x = a\(t), the rate law 
has the form dc/dt = —k(t)(a — x) where kit) is the instantaneous rate constant and 
this form implies 

If Aoo is known from accurate experiments or from our computed estimates, then kit) 
is determined; the variation of kit) provides crucial information concerning reaction 
kinetic mechanism and energetics, from at least one theory recently developed for ele- 
mentary reactions [3 1 and for such theories and developments, it may be anticipated that 
nonlinear methods would be used to accurately determine kit) that would yield the so- 
called "reactivity coefficients" [3 | that account for variations in k that would provide 
fundamental information concerning activation and free energy changes. Figure© 
refers to the computations under the assumption of first order linearity of concentration 
and the conductivity. Whilst very preliminary, non-constancy of the rate constants are 
evident, and one can therefore expect that another area of fruitful experimental and 
theoretical development can be expected from these results. 



2 2. CONCLUSIONS 

The results presented here provides alternative developments based on NLA that is able 
to probe into the finer details of kinetic phenomena than what the standard represen- 
tations allow for, especially in the the areas of changes of the rate constant with the 
reaction environment. Such studies would involve building up another set of axioms 
that is consistent with a varying k(t) kinetic coefficient. Even with the assumption 
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Figure 5: Variation of k with time or concentration changes based on the experimental 
value Aqo = 2050units and the computations based on different polynomial degrees 
n = 2, 3, 4 and the computed Aoo values for Method 1 and Method 2 for fixed polyno- 
mial degree n = 3. 

of invariance of k, one can always choose the best type of polynomial order that is 
consistent with the assumption, and it appears that the initial concentration as well as 
the rate constant seems be be predicted as global properties based on the polynomial 
expansion. 
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